Broadband chaos generated by an opto-electronic oscillator 
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We study an opto-electronic time-delay oscillator that displays high-speed chaotic behavior with 
a flat, broad power spectrum. The chaotic state coexists with a linearly- stable fixed point, which, 
when subjected to a finite-amplitude perturbation, loses stability initially via a periodic train of 
ultrafast pulses. We derive approximate mappings that do an excellent job of capturing the observed 
instability. The oscillator provides a simple device for fundamental studies of time-delay dynamical 
systems and can be used as a building block for ultra-wide-band sensor networks. 
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A deterministically chaotic system displays extreme 
sensitivity to initial conditions and the spectra of the fluc- 
tuating system variables are broadband. Yet, for typical 
chaotic devices, the power spectra often contain several 
sharp features that stand out above a broad background, 
which are often associated with weakly unstable peri- 
odic orbits that are part of the backbone of the strange 
attractor. The fact that the power spectra for typical 
chaotic devices are not featureless limits their applica- 
tion in ultra-wide-band (UWB) sensor networks 1] and 
in chaos-based ranging devices [2], for example. 

In this Letter, we describe an opto-electronic time- 
delay oscillator that displays high-speed chaos with an 
essentially featureless power spectrum. The chaotic be- 
havior coexists with a linearly stable quiescent state. If 
the system starts in this state, a finite-size perturbation 
is needed to force the system to the chaotic state. We 
show that a sufficiently large perturbation causes the sys- 
tem to produce an initially periodic train of ultrafast 
pulses whose spacing and amplitude becomes irregular 
for longer times. Our observations are in good agreement 
with the predictions of a nonlinear stability analysis of 
the fixed point. 

Our work has important implications for understand- 
ing the stability of general time-delay systems, for which 
coexisting states are common. For example, the stability 
and noise sensitivity of opto-electronic microwave oscil- 
lators [3], synchronized neuronal networks [4(, synthetic 
gene networks and controlled chaotic systems [6, 7] 
may be adversely affected by the presence of a coexist- 
ing chaotic state. Our analysis predicts the amplitude of 
noise or externally applied perturbations that allow such 
systems to 'sense' the coexisting strange attractor. 

Our opto-electronic oscillator consists of a nonlinear el- 
ement placed in a time-delay feedback loop and displays 
a variety of dynamical behaviors that depend on system 
parameters. As shown in Fig.[TJ the beam generated by a 



semiconductor laser (wavelength 1.55 /im) is injected into 
a single-mode optical fiber, passes through a polarization 
controller, and is injected into a Mach-Zehnder modula- 
tor (MZM). The transmission of the MZM is a nonlinear 
function of the applied voltage, where we independently 
apply a time-dependent voltage to the radio-frequency 
(RF) port of the device (half- wave voltage V 7V} rf=7A V) 
and a dc- voltage Vb to bias it at any point on the trans- 
mission curve (half- wave voltage K-,dc=7.7 V). Light ex- 
iting the modulator passes through an additional piece of 
single-mode fiber (length ^5 m) serving as a delay line 
and is incident on a photodetector. Half of the result- 
ing signal, denoted by V, is amplified by an inverting 
modulator driver (gain gMD = —22.6) and fed back to 
the MZM via the ac-coupled input port. The modulator 
driver saturates at high voltage with saturation voltage 
V sa t = 9.7 V. The other half of the signal is directed to 
a high-speed oscilloscope (8 GHz analog bandwidth, 40 
GS/s sampling rate). The total delay of the feedback 
loop T=24.1 ns. 
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FIG. 1: Experimental setup. Inset: Nonlinear transmission of 
the MZM (ratio of the output to input powers of the device) 
as a function of the dimensionless bias voltage m. 

Similar opto-electronic oscillators have been studied 
previously, dating back to the seminal work of Ikeda @] • 
One distinguishing feature of our device is that the am- 
plifier is ac-coupled so that feedback of low frequencies 
is suppressed. Also, high frequencies are suppressed due 
to the finite response time of the photodetector and am- 
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plifier. We find that the linear frequency response of the 
various components of the system is well described by 
a two-pole bandpass filter with a low- (high-) frequency 
cut-off uj- = 1.5 x 10 5 s" 1 (uj + = 7.5 x 10 10 s" 1 ), center 
frequency loq = ^/uj_Ul! + = 1.1 x 10 8 s _1 , and bandwidth 
A = uo+ — uj- = 7.5 x 10 10 s -1 . The system is described 
in terms of a single time-delay integro- differential equa- 
tion (or, equivalently, two coupled time-delay differen- 
tial equations) as opposed to a single time-delay differ- 
ential equation used by Ikeda. Such time- delay integro- 
differential equations display unique bifurcation proper- 
ties 0, LL0| and new behaviors, such as chaotic breathers 

Another important distinction of our work is that we 
bias the MZM at the maximum of the transmission curve 
(m = ttVb/^Vtt^c = 0) shown in the inset of Fig. [U 
As shown below, such a bias renders the quiescent state 
of the system linearly stable. Essentially all other re- 
search has focused on the case where the bias is set 
to the half-transmission point of the transmission curve 
(m = 7r/4), where the quiescent state is most linearly 
unstable. Counter examples exist, such as the work of 
Meucci et al. [12|, although they did not focus on the 
behavior discussed here. 

The dynamics of our opto-electronic oscillator is de- 
scribed by the dimensionless coupled time- delay differen- 
tial equations (DDEs) [To|, [Tl| 

x(s) = —x(s) — y(s) + c[x(s — t)] (1) 
y(s) = ex(s). (2) 

Here, c[x] = 7 cos 2 (m + d tanh x) —7 cos 2 m is the nonlin- 
ear delayed- feedback term, x = gMDV/V sa t, the overdot 
denotes the derivative with respect to the dimensionless 
time s = £A, 7 is the overall gain of the feedback loop 
and is proportional to the laser power, r = TA, and 
e = uJq/A 2 characterizes the bandpass filter. Differing 
from Ref. [llj, we incude a hyperbolic tangent function 
in c[x] to account for amplifier saturation, characterized 
by the parameter d = TtV S at / IV-* ,rf • In our experiments, 
three parameters are held fixed (d = 2.1, r = 1820, and 
e = 2.0 x 10~ 6 ), while 7 can range from — 5 by adjusting 
the laser power with an attenuator and m ranges from 
— 7r/2 to 7r/2. For future reference, note that x and V 
have opposite signs because gMD < 0. 

We first investigate the linear stability of one of the 
fixed points of Eqs. (pQ) and (0 located at (x*,y*) = 
(0, 0), which is the quiescent state of the oscillator. Lin- 
ear stability analysis predicts that the fixed point is stable 
for small 7 and undergoes a Hopf bifurcation (a transition 
to an oscillatory behavior) at 

1H = ~d^2mj' (3) 

where b± is a constant that depends on r and e and is 
approximately equal to ±1 for our conditions. Clearly, 
the fixed point is linearly-stable for all 7 for m = where 
7h diverges. 



Experimentally, we find that the stability properties of 
the fixed point are much more complex than predicted by 
the linear analysis presented above. In particular, near 
m = 0, we find that finite-size perturbations destabilize 
the fixed-point for 7 < which can only be understood 
from a global (nonlinear) stability analysis of the model. 

In the experiments, we slowly increase 7 from zero until 
the fixed point (x*,y*) = (0,0) is destabilized. It is seen 
in Fig. [2^ that there is very good agreement between the 
predictions of the linear theory and experiment around 
m = ±7r/4 (the standard bias used in most previous ex- 
periments), but there is substantial disagreement in the 
vicinity of m = (see Fig. [2Jd). At m = 0, the system 
loses stability by transitioning directly to a broadband 
chaotic state at 7 = 4.36, undergoing transient pulsations 
en route to chaos. The finite-size perturbations needed 
to destabilize the fixed point originate from noise in our 
system (e.g., laser relaxation oscillations, and detector 
dark and shot noise). 
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FIG. 2: Observed values of 7 for which the system transitions 
from steady-state to oscillatory or pulsing behavior as a func- 
tion of m, with 7_ff superimposed (solid line). The squares 
in (a) and (b) indicate low experimental noise, while the dia- 
monds in (b) indicate a higher level of noise. 

We add more noise to the system to measure its af- 
fect on the instability threshold. Insterting into the loop 
an erbium-doped fiber amplifier (EDFA) followed by an 
attenuator so that the total optical power injected into 
the oscillator is the same, we observe that the fluctuating 
part of the voltage V increases by a factor of 2.3 (root- 
mean-square noise over a bandwidth from dc to 8 GHz). 
The open diamonds in Fig. [2] show that the instability 
threshold decreases for \m\ < 0.1 due to the increased 
noise. There is also a pronounced asymmetry in the in- 
stability threshold about m = 0. 

In the vicinity of m = 0, we observe that the sys- 
tem loses stability by generating a sequence of ultrashort 
pulses spaced initially by T with a pulse duration (full 
width at half maximum) of ~ 200 ps. Typical transient 
behavior is shown in Fig. |3K, where we have removed the 
EDFA from the setup and adjust 7 just above the insta- 
bility threshold. To more carefully study this transient 
behavior, we add an additional 3-dB power splitter to the 
feedback loop, lower 7 so that the system is in the qui- 
escent state and inject 200-ps-long electrical pulses into 
the loop. For small pulse amplitude, the perturbation 
decays. For sufficiently large pulse amplitude, we ob- 
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serve that the perturbation grows rapidly initially, levels 
off, and the waveform becomes more complex, similar to 
that shown in Fig. [3^. The open triangles shown in Fig. 
|3t> give the critical value of the pulse amplitude needed 
to destabilize the fixed point as a function of the feedback 
loop gain. 
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FIG. 3: (a) Experimentally observed transient behavior that 
results for m = and 7 — 4.36 when the system leaves the 
steady-state. The pulses have a FWHM ~ 0.2 ns and are 
separated by the time-delay T. (b) The critical pulse am- 
plitude as a function of 7 in the experiment (triangles) and 
simulation (stars) with the unstable fixed point of the map 
superimposed. 

The instability boundaries shown in Fig. [2] and the 
transient pulsing behavior shown in Fig. [3] can be under- 
stood by asymptotic analysis of Eqs. (pQ) and ([2]). Con- 
sider a finite-size perturbation to x in the form of a short 
pulse of amplitude — x centered at time s = 0. Due to 
the smallness of e ~ 10~ 6 , y responds extremely slowly 
to this pertubation, and hence we set y = 0. The pertur- 
bation propagates through the feedback loop and, near 
time r, the feedback term begins to grow from zero, cor- 
responding to the presence of the pulse a time r earlier. 
Approximately, this produces another pulse c[xq] in x. 
Thus, this second pulse at time r will generate a third 
pulse at time 2r, and so on, consistent with the behavior 
shown in Fig. [3^. 

From the asymptotic analysis, we derive a one- 
dimensional discrete map of the form 



(4) 



where Xyi can be thought of as the amplitude of a pulse 
at time nr. Further analysis reveals that there exists a 
continuous mapping of a time interal of length r into the 
next future interval of length r. It is given by 



(5) 



where x a ^ denotes the temporal evolution of x over the 
interval (a, b]. One should keep in mind, however, that 
these mappings only approximately predict the dynamics 
of the physical system, as reducing the coupled DDEs to 
a mapping erases all of the effects of the bandpass filter. 

The discrete map ([4j) has multiple fixed points, depend- 
ing on the values of 7, m, and d. Here, we focus on the 
case of m = 0, where there are one or three fixed points 
for d = 2.1. We find that the fixed point x* x = is al- 
ways stable. It corresponds to the steady state where no 



pulses are generated. The other two fixed points emerge 
at a critical feedback gain j c = 0.73, exist for 7 > 7 C , and 
are both negative. The fixed point with the smaller mag- 
nitude, denoted by X ^ . IS unstable, while the fixed point 
with the greater magnitude, denoted by x* 2 , is stable. It 
corresponds to a periodic pulsating state with amplitude 
x* 2 . Thus, the critical perturbation size is given by |x* | 
because perturbation amplitudes greater than this value 
will grow in time towards the stable fixed point x* 2 . For 
7 > 1, where |x*| <C 1, a very good approximation is 
given by x* u « -l/^d 2 ). 

We determine x* numerically from map ([Hand convert 
to physical units (see Fig. [3b)- We also determine x* from 
a numerical simulation of Eqs. ([1} and ([2]). It is seen 
that the agreement between the experiments, and pre- 
dictions of the mapping and simulations are very good. 
Most importantly, it is seen that the minimum pertur- 
bation size decreased as a function of 7, implying that 
noise will eventually destabilize the fixed point for suf- 
ficiently large 7. A similar procedure can be used to 
determine the threshold gain 7^ required to destabi- 
lize the fixed point x* sl = for a given noise intensity 
D = y/2 < x 2 >. At threshold, < x 2 >= |x*| 2 , which 
yields j th ~ l/[dsm(Dd/y/2 - 2m)]. 

For the whole range of m, 7 C continues to indicate the 
birth of two fixed points, which we determine numerically 
from map ([4]) and display in Fig. [4^. Pulsing behavior 
is possible for 7 > 7 C . There is a strong asymmetry in 
7 C about m = 0, indicating that pulsing behavior is least 
likely around m = 7r/4. Also shown in Fig. H^t are 7^ 
and jth for one value of the noise intensity. For 7^ < ^ H 
(~ — 7r/4 < m < 0.1), we predict that the fixed point 
x* sl = tends to be destabilized by a pulsing instabil- 
ity. For < 1th, 7 ~ and small noise, we predict 
that the fixed point tends to be destabilized by the Hopf 
bifurcation determined from the linear stability analysis. 
Thus, we predict that the fixed point will be unstable 
for 7 > min[7#, 7^]. We see that there is qualitative 
agreement between min^, 7^] shown in Fig. HJi and 
the high- noise experimental measurements (Fig. [2Jd). 
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FIG. 4: (a) Instability thresholds using D = 0.28 in j th . (b) 
Temporal evolution of an initial Gaussian pulse (solid line) 
after one iteration of the continuous mapping for m = with 
7 = 1.23 and xo = —0.2 (inner dashed line), and for m = 7r/4 
with 7 = 0.71 and xo = —0.35 (outer dashed line). 

Analysis of the map only gives information about the 
pulse peak amplitude; information about changes in the 
pulse shape is predicted by the continuous mapping ([5]). 
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FIG. 5: The experimental time series (a) and power spectral 
density (b) of the broadband chaotic behavior in the physical 
system for m = and 7 = 4.80 (upper trace). The power 
spectral density of the noise floor obtained for m = and 
7 = 4.30 (lower trace) is also shown. Numerical time series 
(c) and power spectral density (d) for m = and 7 = 4.80. 



In particular, we take the intial perturbation to be a 
Gaussian pulse centered at s = 0, shown in Fig. |4)d. 
For m = 0, we find that the pulse remains symmetric 
and undergoes substantial pulse compression after one 
round trip through the loop (inner dashed curve). This 
pulse compression continues each iterate through the 
loop, eventually, becoming so short that our asymptotic 
analysis breaks down when the pulse spectrum exceeds 
the bandpass filter width. Thus, we predict that per- 
turbations tend to produce ultrafast pulses whose spec- 
trum fills the available device bandwidth. Strong fre- 
quency mixing, known to occur in time-delay systems, 
then gives rise to pulse-to-pulse coupling and chaos [13| . 
For m = 7r/4, the pulse also remains symmetric, but the 
pulse width expands in a round- trip through the loop. 
Thus, short-pulse perturbations tend to smooth out and 
the system will tend to display more sinusoidal behav- 
ior near the instability threshold. For other values of 



771, the pulses become asymmetric and pulse compres- 
sion (expansion) occurs when 7 C < 7^ (7^ < j c .) In 
the experiment, we observe that, in general, transient 
behavior similar to that in Fig. [3^ eventually gives way 
to chaotic behavior such as the one shown in Fig. [5^ for 
the case m = and 7 = 4.80, which is just above the 
noise-induced threshold seen in Fig. El Figure [5Jd shows 
the one-sided power spectral density (PSD) of the exper- 
imental chaotic time series with a resolution bandwidth 
of 8 MHz. Interestingly, the power spectrum is essen- 
tially 'featureless' - flat up to the cutoff frequency of the 
oscilloscope used to measure the dynamics (8 GHz). In 
greater detail, the spectrum is contained with a range 
of 15 dB with a standard deviation of 3 dB for frequen- 
cies below 8 GHz. Our observation is consistent with 
our analysis of the continuous map above, and indicates 
that all frequencies contribute with approximately equal 
strength and that there are no weakly unstable periodic 
orbits embedded in the strange attractor. This behavior 
contrasts with most other chaotic systems where numer- 
ous large peaks appear in the power spectrum. 

We compare our results to the case where 7 = 4.30, 
which is just below the threshold for noise-induced insta- 
bility. As seen in Fig. [5Jd, the power spectral density is 
at least 40 dB below that observed when the oscillator 
is in the broadband chaotic state and is consistent with 
the noise-floor of the overall system. Furthermore, the 
noise floor is contained within a range of 18 dB with a 
standard deviation of 2 dB, indicating that the chaotic 
spectrum is nearly as featureless as the spectrum of the 
system noise. Also, once the system exhibits broadband 
chaos at m = for sufficiently high gain, the dynamics 
are not sensitive to changing the bias point by 1 V either 
way. Figures |5l:,d show similar broadband behavior of 
the numerical solution of Eqs. (1) and (2). 
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